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ATTRACT 

This  report  presents  a  simplified,  semi- theoretical  model  of 
fractionation,  suitable  for  making  interim  estimates  of  degree  of 
fractionation  and  of  radionuclide  partition  between  local,  inter¬ 
mediate,  and  worldwide  fallout.  The  principles  set  forth  are  applicable 
to  the  treatment  of  air-,  tower-,  and  surface-burst  debris  (in  the 
order  of  decreasing  confidence)  and  to  correcting  fallout-prediction 
systems  for  fractionation  effects.  The  material  provides  the  first 
step  necessary  to  illustrate  theoretically  the  definition  of  contami¬ 
nation  level  proposied  in  Part  II  of  this  series. 
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SUMMARY 


Problem 

The  effects  of  radionuclide  fractionation  severely  complicate 
the  prediction  of  many  properties  of  nuclear' tomb  debris. 

Findings 

A  semi- theoretical  model-  can  be  used  in  a  relatively  simple 
manner  to  estimate  both  radionuclide  fractionation  and  its  effects 
on  exposxire  dose  rate.  The  model  is  recommended  for  illustrative 
purposes,  rule  of  thumb  estimates,  and  as  a  stop-gap  until  either 
better  models  or  m.ore  extensive  information  becomes  available. 
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PREFACE 


This  series  of  reports  presents  and  discusses  the  effects  of 
radionuclide  fractionation,  in  nuclear  bomh  debris.  Part  I  (Ref.  3) 
defined  fractionation  as  "any  alteration  of  radionuclide  composition 
occurring  between  the  time  of  detonation  and  the  time  of  radiochemical 
analy  is  which  causes  the  debris  sample  to  be  nonrepresentative  of 
the  detonation  products  taken  as  a  whole."  It  showed  how  the  radio¬ 
nuclide  compositions  of  fractionated  samples  could  be  correlated 
empirically  by  logarithmic  relations.  Part  II  (Ref.  l)  used  these 
relations  as  the  basis  of  a  technical  discussion  of  contamination 
density  as  applied  to  fractionated  nuclear  debris.  This  part  presents 
a  theoretical  foundation  for  the  observed  logarithmic  correlations 
of  Part  I.  It  uses  this  as  a  ^Impliried  of  c'timati'ig  frac¬ 

tionation  a.s  a  function  of  particle  size  and  the  partition  of  product 
radionuclides  among  local,  intermediate,  and  worldwide  fallout.  As 
an  interim  prediction  system  it  is  applicable,  with  decreasing  con¬ 
fidence,  to  air- burst,  tower-burst,  and  surface-burst  debris.  Part  IV 
will  extend  the  calculations  to  show  how  fractionation-correlation 
parameters  can  be  used  to  estimate  the  exposure  dose  rate  from  nuclear 
debris  with  various  degrees  of  fractionation.  It  will  serve  to 
illustrate  the  proposals  made  in  Part  II. 
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INTRODUCTION 


Part  II  in  this  series^  recomnended  a  definition  of  surface 
density  of  fallout  contamination  which  was  realistically  related  to 
the  spatially  variable  radiochemical  composition  in  fallout  patterns. 
In  order  to  illustrate  the  principles  presented  in  that  report,  one 
needs  a  means  of  predicting  the  quantities  to  be  used  in  the  illus¬ 
tration.  One  means  of  mjdting  such  predictions  is  found  in  a  model 
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proposed  by  C.  F.  Miller.  However,  Miller's  model  is  very  complex, 
it  presents  conceptvial  difficulties,  it  requires  machine  computation 
for  its  employiaent,  it  needs  input  data  which  is  presently  unknown, 
and  it  is  in  need  of  considerable,  fundamental  revision  (cf.  Ref.ll). 
There  is  a  definite  need  for  a  simplified,  interim  model  for  pre¬ 
dicting  the  effects  of  radionuclide  fraoti.onation  on  radiation  hazard 
from  fallout. 

This  report  describes  a  conceptually  simple,  easily  usable,  semi- 
theoretical  model  which  is  applicable  to  air,  tower,  and  land- surface 
burst  debris.  It  shows  the  relation  of  the  model  to  observed  cor¬ 
relations  of  fractionated  debris  and  uses  this  relation  to  predict 
radionuclide  partition  between  local,  intermediate,  and  worldwide 
fallout. 
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In  Part  IV  of  this  series,  this  model  will  fona  a  "basis  for 
illustrating  calculations  of  exposure-dose  rate  according  to  the 
newly  defined  surface  density  of  contamination. 

1 

Familiarity  with  the  preceding  reports  of  this  series  ^  is  pre¬ 
requisite  to  understanding  the  material  presented  here. 


THEORY 

The  mtdel  to  "be  developed  rests  upon  two  principal  assumptions. 
The  first  is  that  the  nuclear  dehris  consists  of  tnacroscopically 
homogeneous,  spherical  particles  with  a  lognormal  size- frequency  dis¬ 
tribution.  The  second  is  that  the  ultimate  distri'bution  of  each  mass 
chain  anong  the  particles  is  proportional  to  some  power  of  the  par¬ 
ticle  diameter.  Additional  assumptions  will  he  required  to  apply 
the  model  to  typical  situations. 

In  this  section  we  will  first  aescrihe  ■ciic  iognori;::;..  dtstra’vution 
and  its  pertinent  features.  We  will  then  apply  these  features  to 
the  interrelation  of  size,  surface,  and  volume  for  a  lognormal  distri¬ 
bution  of  spherical  particles.  Finally,  the  relation  of  these  features 
to  fractionation  is  introduced  by  means  of  the  second  assumption.  A 
relationship  "between  pairs  of  mass-chain  ratios  is  thereby  o"btained 
which  is  independent  of  particle  size  and  in  agreement  with  observed 
properties  of  fractionated  debris. 
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The  Lognormal  Distribution 

The  l.ognormal  distrihution  is  frequently  used  to  describe  nuclear 

4 

debris.  Stewart  derived  a  lognormal  particle  size- frequency  distri¬ 
bution  from  theoretical  considerations,  but  his  derivation  is  much  less 
realistic  when  applied  to  a  land- surface  burst  than  it  is  for  an  air 

burst  or  tower  burst.  Anderson^  has  adopted  a  lognormal  distribution 

* 

of  radioactivity  with  particle  size  for  use  in  bis  D-Model.  Miller 
has  proposed  a  lognormal  distribution  of  mass  with  particle  size  for 
use  in  his  fractionation  model  for  land-surface  bursts.  It  will  be 
shown  below  that  these  proposals  are  in  haimony.  Although  lognormal 
distributions  of  various  types  are  frequently  observed  in  samples  of 
fallout,  it  remains  to  be  established  whether  these  distributions 
apply  to  the  total  amount  of  debris  produced. 

A  random  variable  is  lognormally  distributed  if  the  logarithm  of 
the  random  variable  is  normally  distributed.  Thus,  x  is  lognormally 
distributed  with  mean  n  and  variance  if  its  frequency  distribution 
function  is 


dA(u,  0^)  _  1  exp 

d  In  X  o'^/St 


1 

2 


In  X  -M 


(1) 


The  notation  is  essentially  that  of  Aitcheson  and  Brown.  The  quantity 
d  A (p,  o^)/d  in  X,  when  multiplied  by  d  In  x,  gives  the  probe-bility 
that  a  randomly  chosen  sample  has  a  value  of  x  whose  logarithm  lies 


*C.  F.  Miller,  Private  Communicaii^u. 
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between  In  x  and  In  x  +  d  In  x.  Appendix  A  of  this  report  gives 

elementary  background  material  on  probability  distributions  and  their 

relation  to  size-frequency  distributions. 

An  interesting  and  useful  property  of  the  lognormal  distribution 

is  the  following  (cf.  Ref.  6,  Theorem  2.6).  The  kth-moment  distribution 

of  a  A  distribution  with  mean  and  variance  is  also  a  A  distribution 
2  2 

but  with  mean  H  +  ko  and  variance  o  .  Thus,  it  is  easily  verified 


that 


d  In  X 


k 

=  -12 


d  ACiig,  0^) 

d  In  X 


(2) 


where  x^^  =  (x^^  x^^  =  exp  A  useful  form  of  Eq,  2  is 

2 

obtained  by  writing  =  n  +  k^o  .  Then 


d  A(ja  +  0^  o^)/d  In  X  Jx 

d  A(n  +  kp  0^,  0^)/d  In  X  |  [,  \  2  | 


l'l  V^l 

J. 


(3) 


The  Interrelation  of  Size,  Sturface,  and  Voiur.'ie  DisLji'.jt V.ns 

As  a  consequence  of  Eq.  2  above,  if  the  frequency,  the  surface, 
or  the  volume  of  a  collection  of  spherical  particles  is  lognormally 
distributed  with  respect  to  the  particle  size,  the  other  two  quantities 
are  also  lognormally  distributed  with  respect  to  the  particle  size.  For 
exemp.le,  let  x  represent  the  particle  ditimeter  and  let  its  normalized 

frequency  distribution  function  be  d  A  (iJi^,  _  If  N  is  the  total 

d  In  X 

number  of  particles  and  S  the  total  surface  of  the  particles,  the 
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distributicn  of  the  surface  with  respect  to  x  is 


d  In  X  d  In  X  ' 

Here  Xg  =  2n  (2o^),  so  =  x^,  exp  (o  ). 

Since  ^  is  noiinalized,  the  total  surface  is 


S  =  It  N  x^  . 

"hs 

Similarly,  the  distribution  of  volume  with  x  is  given  by 

,x3»  d  A(v  O^) 

6  d  In  X  6  d  In  X 

«  y  X  =  X  exp  (30^)  and  the  total  volume  is 
—v  Hi 

V  =  «  N  sL 


(5) 

(6) 


(7a') 


S  X 
=  5  -sv 

Thc  total  mass  of  the  particles  is 

M  =  pV 

where  p  is  the  density.  The  mass  is  obviously  distributed  like 
the  volume  since  the  particles  are  homogeneous,  and-  the  mean  of  the 

mass  distribution  (1^^^  equals  \x^,  and  consequently 

An  important  quantity  for  our  purposes  will  be  the  ratio  of  the 
fraction  of  the  total  surface  to  the  fraction  of  the  total  volume  in 


(7b) 

(8) 


a  given  particle-size  increment  as  a  function  of  x.  From  Eq.  3  above,' 
this  ratio  is 

cL  A(u  ,  o^)/d  In  X 

s,v  5 - 

d  0  )/d  In  X 

=  exp  .(5d^/2) 

,  X  '  / 

“sv/x  (9a) 

=  6  V/Sx  (Qb) 


According  to  Eq.  9a,  r  is  unity  when  x  is  equal  to  the  geometric 

6,V 

mean  of  x  and  x,  . 
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RELATION  OP  THE  LOCWORMAL  DISTRIBUTION  TO  FRACTIONATION 


Consider  now  two  mass  chains  i  and  j  which  are  distributed  among 

k,  k, 

particles  of  fallout  of  diameter  x  according  to  x  and  x  respec¬ 
tively,  the  fallout  particles  uhemselvec  naVi.^^  friq  "-"ics  vhlch  are 
lognormal  in  x.  The  ratio  of  the  fraction  of  mass  chain  i  to  the 
fraction  of  mass  chain  j  in  particles  of  size  x  is  found  from  Eq.  3 
to  be 


d  +  k^  a^,  o^)/d  In  x 

d  +  k,  0^  a^)/d  In  X 

n  V 
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'I’wo  ratios  for  two  piiirs  of  mass  cKains  can  be  related  by  eliminating 
Xj^/x  from  their  respective  equations: 


l/(kj  ^  k^) 


I  ^  ■ 


l/(k  -  -k  ) 
'  V  u' 


exp  I  - 

L  ■ 


k  +  k 
u  V 


k.  -  k. 


k  -  k 
V  u 


Vv  *  (‘'J  -  “l)  \ 


This  equation  indicates  a  linear  relationship  between  the  logarithms  of 

mass-chain  ratios  which  is  independent  of  particle  size.  Linear  relation¬ 
's 

ships  between  logarithms  of  mass-chain  ratios  have  in  fact  been  observed. 


APPLICATION  OF  MOliEL  TO  LAND  SURFACE  BURSTS 


We  will  now  illustrate  the<  foregoing  model  by  applying  it  to  some 
land-surface  fallout  distributions  currently  in  use  at  NliDL.  We  will 
compare  these,  first  according),  to  their  physical  properties,  then,  with 
the  help  of  some  additional  assumptions,  according  to  their  predicted 
radiochemical  properties.  We  will  be  particularly  concerned  with  the 
•variation  of  fractionation  •with  particle  size  and  with  the  partition 
of  radionuclides  among  local,  intermediuLe,  and  worldwide  fallout. 

The  Distributions 

Three  lognormal  distributions  have  been  used  at  NRDL  to  describe 
fallout : 
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The  H-61  distrihution  is  used  for  surface  bursts  on  Nevada  soil 
regardless  of  yield.  According  to  this  distribution,  half  of  the  gross 
activity  in  the  nuclear  clovid  is  in  particles  vith  cianeters  of  100  u  or 
less  and  1  >  is  in  particles  with  dianeters  of  10 '  ti  or  greater. 

The  C-6l  distribution  is  used  for  surface  bursts  on  coral  surfaces. 

It  is  also  yield- invariant.  According  to  this  distribution,  half  of  the 
gross  activity  is  in  particles  with  dianeters  of  100  n  or  less  and  i  ^  is 
in  particles  with  dianeters  of  5000  (i  or  greater. 

Eie  distribution  proposed  by  Miller^  contains  half  of  the  rriass 
in  particles  with  diaseters  of  100  51  or  less  and  0.4  ^  of  the  aiass 
in  particles  with  dianeters  of  10^  or  nore  (l  ^  in  articles  of  diameter 
5600  n  or  greater). 

In  order  to  illustrate  the  dependence  of  various  distribution 
properties  pn  the  mean  and  variance,  ve  will  treat  the  N-61  and  C-61 
distributions  as  thovigh  they  applied  to  mass  (or  vclvce,  since  the  par¬ 
ticles  are  assumed  homogeneous)  msusad  of  to  sjr&ss  ic  r~.ty.  Crmspsrison 
of  the  three  distributions  on  this  basis  will  illustrate  the  sensitivity 
to  the  choice  of  variance.  To  illustrate  the  sensitivity  to  the  choice 
of  mean,  we  will  use  two  distributions  with  the  same  variance  as  that 
proposed  by  Miller,  but  with  mass-modal  particle  sizes  of  50  snd  200  ji. 

Ne  will  designate  these  distributions  end  respectively,  and  the 

original  MUer  distribution  by 

Table  1  summarizes  some  important  statistical  and  physical  properties 
of  these  distributions,  together  with  the  equations  involved  in  their 
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Table  1 


Properties  of  Various  Distributions,  Illustrating  Their 
Sensitivity  to  Metui  and  Variance 


Quantity 

c-6i 

U-61 

«100 

^■^200 

Variance  (Dimensionless) 

1.682 

1.980 

1.736 

1.736 

1.736 

a 

2.829 

3.919 

3.015 

3.015 

3.015 

Modal  Diameters  (u) 

X 

—V 

100 

100 

50 

100 

200 

X  =  X  exp  (-0  ) 

-s  -V  ^  ' 

x^  =  (-30  )  2 

5.907 

.04x10"^ 

1.986 

7.81x10"^ 

2.452 

5. 87x10“ 2 

4.904 

1.17x10"^ 

9.809 

2.35x10"' 

X  =  ■\l  X  X 

— sv  y  —s  -V 

24,3 

i4.i 

n.i 

22.1 

44.3 

Means  (Dimensionless) 

=  In  X 

V 

4.605 

4.605 

3.912 

4.605 

5.298 

(1  =  In  X 

1.776 

0.686 

0.897 

1.590 

2.283 

(Ji  =  In  X 
n  -n 

Physical  Properties 
S/V  =  6/Xg^  (H”^) 

-3.876 

-7.151^ 

-5.137 

-4.449 

-3.751 

0.247 

0.426 

0.541 

0.271 

0.3.35 

N/V  =  6/«  Xn^3 

0.655 

87.5 

12.0 

1.51 

0,187 

Fractions  of  Mass  or  Volume  in 

Particle-Size 

Ranges 

0-25  ^ 

0.205 

0.242 

0.345 

0.212 

0.115 

25-50  \i 

0.135 

0.121 

0.155 

0.133 

0.097 

>50  \i 

0.660 

0.637 

0.500 

0.655 

0,788 

Fractions  of  Surface 

in  Particle-Size  Ranges 

0-25  \i 

0.805 

0.900 

0.910 

0.826 

0.705 

25-50  H 

0,093 

0.048 

0.049 

0.084 

0.121 

>50  \i 

0.102 

0.052 

0.051 

0.090 

0.174 

Fraction  of  Surface/Fraction  of  Volume  (r 

S  >  V 

)  for  Particle-Size  Ranges 

0-25  H 

3.93 

3.72 

2.64 

3.90 

6.13 

25-50  \X 

0.689 

0.397 

0.316 

0.632 

1.25 

>50  |1 

0.155 

0.082 

0.102 

0.137 

0.221 
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calculation.  Figure  1  shows  the  distribution  in  differential  form 
(frequency  function)  and  Fig.  2  shows  it  in  integral  form  (the  distribution 


function) . 

We  will  now  apply  The  surface-to-volume  ratio  to  the  estimation  of 
fractionation  with  particle  size,  and  will  use  the  fractions  of  surface 
and  volume  in  various  particle-size  groups  to  estimate  the  partition  of 
individual  radionuclides  among  local,  intermediate  and  worldwide  fallout. 
Predicted  Radiochemical  Composition  vs  Particle  Size 

In  order  to  apply  these  considerations  to  the  estimation  of  radio¬ 
chemical  properties  from  a  land-surface  burst,  two  additional  assumptions 
must  now  be  hade  to  obtain  k  values.  The  first  is  that  the  mass-95  chain 
is  distributed  like  the  volume 
chain  is  distributed  like  the  surface 


|kg^  =  3m  second  is  that  the  mass-89 


^9  “ 


=  2 


.  This  allows  roo  nc  to  be 
o9»95 

equated  to  r^  ^  and  calculated  as  a  function  of  particle  size  from  Eq.  9* 

The  values  of  r  integrated  over  size  increments  are  tetmlated  in 
s,v 

Table  1, 


the  form 


Empirical  correlatiohs  of  radionuclide  ratios  have  been  obtained  in 
3 


^°®10  ^^1,89  “  ®i  \  ^°®10  ^95,89 


(12) 


or 

^°®10  ^i,95  ^^89, 95 

where  a^  and  b^^^  are  regression  coefficients. 

Comparing  Eq.  13  with  Eq.  11  in  view  of  the  above  assumptions, 

J  =  V  =  95,  k.  =  k  =  3,  u  =  89,  k  =2,  Equation  11  becomes 
j  V  u 


(13) 


J-og  + 


(3  -  \)'<\  -2)  o’" 
2  X  2.303 


n'Vt  ^  >»o -f*/*;  v»a 


b.  =  k.  -  2 
i  1 


(l4) 


=  (3  -  k^)(k^  -2)  a^/k,6o6 


=  (1  -  "b^)  a^/4.606 


(15) 


The  evElmtion  of  by  the  radiochemical  analysis  of  fallout 


samples  or  its  estimate  by  the  method  described  in  Ref.  S>  allows 
the  radionuclide  con^xDsitlon  to  be  estimated  as  a  function  of  particle 
size.  If  0  i^i  known,  a^  can  also  be  calculated.  For  a  given  distri¬ 
bution,  woT^d  have  a  maximm  value  for  a  mass  chain  with  b^  equal 
to  1/2,  In  such  a  case 


max 


0^/18.  If 


Predicted  Partition  of  Radimmclides  Among  Lc-oul, 

Intermediate  and  Worldwide  Fallout 

The  assumptions  vised  in  the  preceding  paragraph  (viz.,  kg^  = 
also  allow  estimates  of  the  partition  of  various  radionuclides  among 
local,  intermediate,  and  worldwide  fallout.  For  this  purpose  we  first 
define  local  fallout  as  consisting  of  all  particles  with  diameters  of 
56  n  or  more,  Intennediate  fallout  a"  particles  with  diameters  reinging 
from  25  to  50  H,  and  worldwide  fallout  as  particles  with  diameters;  of 
25.'H  or  less.  The  fraction  of  the  jrai.K.-95  chain  in  any  of  these  portions 
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is  then  equal  to  the  fraction  of  the  total  particle  volume  in  that  por¬ 
tion,  as  shown  in  Table  1.  The  same  is  true  for  any  other  quasi¬ 


refractory  mass  chain  {i.e.,  one  for  which  3.). 

The  fraction  of  the  mass-89  chain  in  any  portion  is  equal  to  the 
fraction  of  the  surface  lying  in  that  portion.  Values  for  the  three 
ranges  cited  are  shown  in  Table  1.  The  same  is  true  for  any  other 
nuclide  for  which  b^  =  0. 

The  fraction  of  intermediately  behaving  chains  can  -be  estimated 
by  calculating  their  modal  volume  values  as  in  Eq.  2: 

r 

=  Xy  exp  |^^(b^  -  1)  (16) 

From  what  has  preceded  it  is  clear  that  the  variance  for  each  such 
distribution  will  be  the  same  as  the  variance  for  any  geometric  par¬ 
ticle  property  for  that  distribution.  Curves  can  then  be  plotted  to 
yield  the  fraction  in  any  group  according  to  the  value  of  b^.  This 
has  been  done  for  the  distribution,  and  the  resai.-  1'*  oho  .Ti  in 

Fig.  3*  'I'he  values  of  b^,  for  some  mass  chains,  as  determined  from 

2^9 

Ref.  3,  are  indicated  on  the  curve.  That  of  Pu  is  estimated. 

APPLICATION  OF  MODEL  TO  AIRBURSTS 

As  mentioned  above,  Stewart's  treatment  of  nuclear  debris  forma¬ 
tion  is  more  realistic  for  air  bursts  than  for  land  surface  bursts. 
His  calculations  apply  to  a  device  composed  of  1000  kg  of  Fe  and,  as 
can  be  shown  from  the  data  he  uses,  a  yxeld  of  10  kt.  The  total 
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The  assvimptions  of  the  previous  section,  namely,  that  the  iiia)!s-95 

chain  is  distributed  among  the  particles  in  proportion  to  their  volume 

and  the  mass-89  chain  is  distributed  according  to  surface,  are  also  more 

realistic  when  applied  to  air  b\u:sts  because  the  debris  remains  in  contact 

longer  with  the  condensing  material.  The  fact  that  fractionated  samples 

7 

of  air  bvtrst  debris  correlate  logarithmically  is  also  in  agreement  with  a 
lognormal  distribution,  as  explained  above.  Therefore,  the  present  model 
can  be  used  also  in  estimating  fractionation  for  air  burst  debris.  However, 
certain  precautions  should  be  indicated.  Thus,  Stewart's  theory  was 
developed  for  bursts  of-  low  yield  with  condensation  times  of  the  order  of 
10  sec.  or  less.  No  theory  for  longer  periods  has  been  found.  However, 
two  modifications  can  be  made  to  Stewart's  formula  which  will  make  it  more 
realistic  and  more  applicable  to  bursts  of  higpier  yield  and  lower  mass. 

The  first  modification  concerns  putting  a  lower  limit  on  the  particle 
size.  With  hi^er  yield- to-mass  ratios,  smaller  particles  in  greater 
abundance  are  to  be  expected.  Ob’/iously,  there  is  a  ixuill.  1..  a  dust-free 
atmosphere,  it  would  be  initially  the  molecular  diameter.  Tropospheric 
aerosols  have  a  lower  limit  of  about  8  x  10  p  diameter  due  to  the  coagula- 

g 

tion  of  smaller  particles. °  Similar  mechanisms  can  be  expected  to  eventvially 
establish  a  lower  limit  of  air-burst  debris  particle  size.  In  a  real 
environment,  the  liicit  would  be  approximately  the  modal  size  of  the  dust 
particles.  In  either  case,  we  can  call  this  limit  x^  and  shift  coordinates 
in  such  a  manner  that  the  distribution  ’.aw  does  not  result  in  smaller 
particles. 
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The  second  modification  concerns  the  variation  of  x  with  the 
device's  ratio  of  mass  to  total  yield  M/W.  The  value  of  will 
vary  as  M/W.  Neglecting  the  effect  of  M/W  on  the  other  parameters 
leads  to  a  direct  proportion  between  x^.  and  M/W.  From  Stewart's 
calculations,  with  M/W  dimensionless, 

10  M/W  (microns)  . 

Introducing  this  and  the  shift  in  abscissa  into  Eq.  I7  gives 


lA 


d  In  X 


-fir 


. 

2 

1 _ 

L  \ 

X  -  X 


\2-l 


(19) 


Figure  4  shows  some  differential  graphs  of  this  function  with  x^  =  10^  la 
for  different  mass-to-yield  ratios.  A  thousand-pound  device  weight  is 
assumed.  Integral  curves  are  shown  in  Fig.  5* 

The  validity  of  these  considerations  end  the  development  of  more 
exact  equations  are  subjects  that  require  further  investigation  and 
lie  beyond  the  scope  of  the  present  report. 


APPLICATION  OF  MODEL  TO  TOWER  BURSTS 


Stewart  considers  Eqs,  17  and  I8  applicable  to  tower  bursts  and 
illustrates  the  application  for  a  t-.-'st  of  the  same  yield  in  which 
25  tons  of  steel  tower  are  vaporized.  For  this  situation  he  calculates 
x^  =  0.2  n.  (Our  approximate  equation  would  give  0.25  |i).  Applying 


15 


the  assumptions  set  forth  above  to  Stewart's  tower  hurst  gives  the 

integral  curves  shown  in  Fig.  6.  As  shown  :hy  the  figure,  the  particle 

size  of  unit  ro„  value  has  decreased  from  the  value  of  x  shown 

-s,v 

in  Table  1  to  a  value  of  2.5  According  to  this  treatment,  therefore, 
a  particle  from  a  tower  burst  should  be  more  depleted  in  Sr®^  than  a 
particle  of  the  same  size  from  a  land-surface  burst. 

DISCUSSION 


This  section  attempts  to  point  out,  with  equal  emphasis,  both  the 
defects  and  the  merit  of  the  present  model.  Although  the  model  is  a 
stop-gap,  whose  raison  d'etre  lies  more  in  necessity  than  in  intrinsic 
worth,  it  will  be  shown  that  its  resemblance  to  reality  is  sufficient 
to  justilY  its  use  for  illustrations,  rule  of  thumb  calculations,  and 
estimates  made  in  ohe  absence  of  reliable  information. 

The  principal  assumptions  involved  in  tjas  devuxi^,;...snt  uivd  applica¬ 
tion  have  been  (i)  a  lognormally  distributed  group  of  homogeneous, 
spherical  particles,  (2)  a  distribution  of  masB-95  chain  among  the  par¬ 
ticles  in  proportion  to  the  volume,  and  (3)  a  distribution  of  the  mass-89 
chain  among  the  particles  in  proportion  to  their  surface.  None  of  these 
assumptions  is  new,  and  only  the  flrsb  is  popular.  The  last  two  can  be 
easily  modified  to  conform  to  other  powers  of  the  diameter,  should  this 
appear  desirable.  In  this  section  to  will  discuss  the  validity  of  each 
assumption  in  turn  and  then  compare  the  present  model  with  the  Miller 
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model.  Finally,  the  extension  of  the  model  to  other  Inirst  types 
(water  surface)  will  he  discussed  and  the  unifying  power  of  the  model 
will  he  explained. 

Validity  of  the  Lognormal  'Distribution 

The  case  for  lognormal  particle- size  distributions  has  already  been 
discussed,  but  it  seems  appropriate  to  summarize  it  here,  Stewart's 
prediction  of  lognormal  distributions  for  nuclear  debris  is  most  validly 
applied  to  air  bursts  and  least  to  land- surface  bursts.  The  meager  data 
availab..e  on  air-burst  debris  is  not  valid  below  diameters  of  5  1^,  but 
data  above  this  size  can  'oe  fitted  with  a  lognormal  form  and  reasonable 
parameters  obtained.  Therefore,  the  data  available  support  this  kind 
of  distribution  but  are  not  sufficient  to  substeintiate  it.  In  air  bursts, 
the  particles  are  extremely  heterogeneous  and  in  land-surface  bursts 
they  are  more  so.  The  properties  observed  must  be  considered  as  an 
average,  and  a  theory  that  can  predict  the  average  does  well. 

In  fallout  samples  from  land- surface  buis  kinds  of  distri¬ 

butions  have  been  observed,  but  the  lognormal  activity-particle  size 
distribution  is  as  good  an  approximation  to  the  observations  as  any. 

It  must  be  emphasized  that  these  samples  have  all  been  biased,  not  only 
by  the  sampling  device  itself,  but  also  by  the  fallout  transport  process. 
Thus,  samples  collected  close  to  ground  zero  were  depleted,  in  small 
particles  and  those  collected  at  greater  distances  were  successively 
depleted  in  large  particles.  » No  representative  sample  of  p&’ticles  from 


-a  3aTiti-Bui:face  burst  hac  ever  buen  obtained,  end  probably  none  ever  «il.l 
be.  TTechnipueB  oT  enal'yeing  leilout  data,  noir  being  annlied  nt  this 
l5abQrator;y„  xiry  lead  to  more  xe£ined  entlmateB  of  the  Initial  portiele- 
slze  distribution. 

'.VaHiidity  oT  the  :x^  JiistrJbutiDn  joT  the  {•ioBBSS  Chain  nnd  the  .y^ 
Jistribution  el  •the  'MaEB-lB9  Qhein 

!Ehe  nBBumption  has  been  made  that  the  -inaBB-95  dhaln  iiiBtrjIbutes 
itselff  throughout  the  thUout  norticQeB  aneording  to  the  eiibe  oT  their 
d±Hme±Er‘i.  in  the  ease  oT  oir-iburEt  c^briB,  t^B  Beems  to  be  n.  reoBon- 
Tsbls  aBBumption.  Air-'burst  dehiiiE  haB  been  3'DunQ  to  be  'veiy  betere- 
•geneouB  -with  regard  to  netbor,  BperaiUn  activity,  densda^,,  dee^  charan- 
teristicB,  and  macroBCDpetc  compDEiifciDn,  and  no  rnlntiDnsMp  ^aaiong  these 
iprppertieB  (has  been  .determined.  lEhe  aetivity  id"  individuEn.  •par.ticleB 
cagpairrs  to  '.vary  aphroxiniat&ly  as  the  nuhe  oT  the  diDmeter  2ar  large 
parMaleE,  and  '.varaeB  progreBHively  rmore  like  the  sguare  oT  the  dinmete" 
as  the  ’particle  aize  decreBPes,  behpvbir  '.n  tn  n’*alitativE  agrae- 

rmerit  with  the  behavior  ejected  trom  the  prcBent  Ttiodel. 

vWith  regard  to  land-aurface  burstB,  two  principal  Mnds  oT  tor- 

n  po 

ticleB  have  been  observed.  ’  One  kind  in  .npherical,  iractionoted, 
erf'  ihigh  apccifCic  act±.vity,,  aad  the  activity  in  diB’trfbirted  ■more  or  Icbe 
regularly  throughout  the  -.volume .  The  ether  kind  in  angular,,  Icbb  trae- 
tionated,  tiT  much  lower  npecilllc  activity,  -and  the  actl.v±ty  is  noncen- 
'trated  on  its  surface.  /Othough  the  latter’  is  numsr.’jus,  the  ioroer 
indkeB  the  larger  contribution  to  the  radiologii^ul  field  nnd  containr 

beon  Isevcnthol,  private  comnunicatlon 
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the  larger  fraction  of  the  inaGs-95  cliain.  The  former  kind  appears  to  he 
comprised  of  particles  that  were  completely  melted  and  perhaps  partially 
vaporized.  The  larger  "bilSGC  sii!rial>ly  ifcH  away  from  nucloat* 
cloud  while  the  raass-89  chain  was  predominantly  in  the  uncondensible 
Kr  form.  The  large,  irregular  particles  are  apparently  formed  by  the 
indiscriminate  condensation  of  vapors  on  cold  particle  surfaces  or  by 
the  scavenging  of  small,  condensed  particles.  Fractionation  was  observed 
to  cause  the  quantity  of  the  mass-89  chain  carried  by  a  particle  to  be 
independent  of  the  particle  type,  by  counteracting  the  increased  activity 
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of  regular  particles  with  Sr  depletion. 

A  further  piece  of  information  on  the  validity  of  the  distributions 
chosen  for  the  maGG-95  and  mass-89  chains  has  recently  appeared  in  a 
report  by  J.  H.  Norman  and  W.  E.  Bell.^^  On  the  basis  of  their  studies 
of  Cs  vaporization  from,  and  condensation  on,  silicate  melts  of  varying 
compositions,  these  authors  conclude:  "Those  eLements  that  exhibit 
relvutively  low  vapor  pressures  wiiji  conclenLi  "■  '."ut  p  'rti  cles  at 

high  temperatures  where  rates  of  diffusion  are  relatively  high.... 

It  is  expected  that  these  elements  will  readily  diffuse  into  fallout 
particles  and  will  become  fairly  uniformly  distributed.  On  the  other 
hand,  elements  that  exhibit  relatively  high  vapor  pressures  will  con¬ 
dense  onto  fal lout  particles  at  low  temperatures,  where  rates  of  diffusion 
are  low.  It  is  expected  that  these  elements  will  be  retained  on  or  near 
the  surface  of  the  particle," 


It  therefore  appears  that,  in  the  absence  of  agglomeration,  there 
are  reasonable  distributions  to  expect,  but  that  further  study  is  re¬ 
quired  to  take  agglomeration  effects  into  accovmt. 

Comparison  vith  the  Miller  Model 


The  Miller  model  and  the  present  one  have  certain  similarities. 

Both  adopt  a  lognormal  mass  vs.  particle- size  distribution  and,  if  the 
particles  are  considered  to  be  completely  melted,  both  will  predict 
that  refiactorily  behaving  mass- chains  are  distributed  among  the  par¬ 
ticles  in  proportion  to  the  cube  of  the  diameter. 

The  Miller  model  has  the  advantages  of  accdmting  for  particles 
which  leave  the  cloud  before  the  ma6s-89  chain  has  condensed  and  of 
being  able  to  handle  partially  melted  particles.  This  is  good  with 
regard  to  the  prediction  of  refractorily  behaving  chains,  but  bad  with 
regard  to  the  prediction  of  volatilely  behaving  chains.  Thus,  the 
closest  that  one  comes  to  accounting  I'or  the  mass-&9  cuai.i  con-ributed 

by  unfractionated,  irregular  particles  is  to  approximate  them  with 
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partialJ.y  melted  spherical  particles  depleted  in  Sr  . 

The  present  model  tends  to  correct  this  defect.  Whether  the  extent 

of  correction  is  so  xmrealistically  large  that  it  is  better  omitted 

remains  to  be  seen. 

Unitive  Power  of  the  Present  Model 

The  applicability  of  the  present  lii'^del  to  air,  tower,  and  land- 

12 

surface  bursts  is  a  point  in  its  favor.  It  lias  been  pointed  out  that 
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the  condensation  processes  occurring  vhen  a  nuclear  device  is  detonated 
on  the  sxu-face  of  a  body  of  fresh  water  can  be  approximated  by  those  of 
an  air  burst  occurring  in  the  presence  of  vaporized  water.  T^.e  present 
model  should  provide  orientation  for  considering  this  type  of  situation. 

In  the  case  of  an  ocean- surface  burst,  the  situation  is  similar  but 
complicated  by  the  presence  of  vaporized  sea  salts. 

All  these  situations  will  receive  consideration  in  due  course. 
Application  to  Models  for  Fallout  Transport 

There  exist  many  models  for  calculating  the  transport  of  radio¬ 
active  fallout  from  the  nuclear  cloud  to  the  ground.  Many  of  these 
consider  the  debris  to  be  divided  into  particles  of  different  sizes. 

By  means  of  the  model  described  here,  the  degree  of  fractionation  can 
be  estimated  for  any  particle-size  range.  The  next  repoi’t  of  this 
series  will  illustrate  the  estimation  of  exposure  dose-rate  as  a  func¬ 
tion  of  time  and  degree  of  fractionation.  With  calculations  of  this 
type,  the  contribution  to  the  dose-rate  from  each  pa- c  Lass  can 
be  estimated  and  fractionation  accounted  foi . 

Resemblance  of  Model  Predictions  to  Observations  of  Nuclear  Debris 

Without  involving  classified  material,  the  resemblance  of  model 
prediction  to  fallout  observations  may  be  summarized  as  follows: 

(1)  The  model  is  in  agreement  with  observed  correlations  of  radio¬ 
nuclide  composition  for  fractionated  debris  in  both  air  and  surface  bursts. 

(2)  Tlie  model  gives  the  observed  trend,  in  radionuclide  composition 
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with  pai-ticle  size  (greater  enrichment  in  volatile-hehaving  chains  with 
decreasing  particle  diameter)  in  both  air  and  surface  hursts. 

(3)  In  a  2-kt,  near-s\irface  hm*st,  the  particle  diameter  for  which 
^89  95  ~  ^  observed  to  be  27  +  17  ji,  in  good  agreement  with  the 
values  of  ^  shown  in  Table  1.  However,  particles  greater  tlian 
several  hundred  microns  were  found  to  be  3.C'SS  depleted  in  volatilely 
behaving  chains  than  the  model  predicts.  This  departure  may  be  due 

to  the  incorporation  of  smaller,  more  volatile-rich  particles  on  these 
larger  particles  by  impaction.  Evidence  of  impaction  also  appeared 
in  microscopic  and  autoradiographic  examinations  of  this  debris. 

(4)  The  model  predicts  reasonable  partition  between  local  and 
worldwide  fallout. 

(5)  The  model  offers  a  significant  improvement  over  the  use  of 
unfractionated  radionuclide  con^xjsitions. 
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PARTICLE  SIZE  in) 


Fig.  1  Freq,UGncy:  Fimction  (Differential)  Curves  of  Mass,  Siurface  and 
Number  of  Particles  as  a  Function  of  Particle  Diameter  for  the  Mj_qq 
Distribution,  The  figure  shows  the  shift  in  modal  value  with  distribu¬ 
tion  moment  and  the  occurrence  of  where  the  surface  and  volume 
curves  intersect. 


NDICfiT-eO 


PARTICLE  DIAMETER  (, 


FRACTION  OF  TOTAL 


WORLD¬ 

WIDE 

fallout 


LOCAL 

fallout 


‘6 


led; 

.00 


FRACTION  OF  TOTAL 


0.01  O.I 


I  10  JO  90  99  99.9  99.99 

PERCENT  OF  PARTICLES  SMALLER'  THAN  SIZE  INDICATED 

Fig.  5  Pariicle  Size  Distribution  -retion  (integral)  Curves  for 
Air  Bursts  of  Various  Mass-to-YleiO  Patios-^ 

*Solld  lines  ore  plotted  for  x-Xq  Dashed  lines  are  plotted  for  : 
iGcumlng  =  10  micron 
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PERCENT  OF  QUANTITY  IN  PARTICLES  LARGER  THAN  SIZE  INDICATED 


PERCENT  OF  QUANTITY  IN  PARTICLES  SMALLER  THAN  SIZE  INDICATED 


Pig.  6  Pai'ticle  3ize  Distribution  Function  (integral)  Curves  for  a 
lO-kt  Tower  Burst  Vaporizing  25  t  of  Steel 
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PARTICLE  DIAMETER  (/l) 


APPENDIX  A 


RELATION  OP  SIZE- FREQUENCY  DISTRIBUTIONS  TO 
PROBABILITY  DISTRIBUTIONS 


If  w^.  let  Nla,b)  be  the  number  of  particles  in  a  given  sample  whose 
diameters  x  satisfy  the  inequality  a  <  x  ^  b,  then  N(0,oc)  is  the  total 
number  of  particles  in  the  sample  and  N(a,b)/N{0,cxo)  =  P{a,b),  the- 
probabilitji  that  a  randomly  chosen  particle  has  a  diameter  x  such  that 
a  <  X  b.  The  size  frequency  function  p(x)  is  related  to  this  proba¬ 


bility  by 


P(a,b)  =  p(x)dx 


J 


(Al) 


If  it  is  desired  to  approximate  the  size- frequency  distribution  with  a 
normal  distribution,  thih.  is  <:.isi.ly  accv-mr.l-'.shed  the  standard  normal 
distribution  function 


y 

by  making  the  equality 


-.4-  exp  (-t^/2) 


from  which 


p(x)dx  =  y  '.It  =  y  dt  dx 
dx 


p(x)  =  y  dt 
ux 


(A2) 


(A3) 


3l 


Thus^  for  a  normal  distribution; 


^  _  1 
dx  ~  a  ^ 


y  == 


exp 


P  (x)  =  .  exp 

ay  2k 


For  a  lognormal  distribution 
.  In  x-^i 


dX  _  1 

dx  ”  xo  ' 


(A4) 


y  = 


exp 


1  1 

f  In  x-p\ 

■  2 

1  "  )  J 

P  (x)  = 


1  1 

In  x-n| 

’  2  1 

a 

=  (A5) 

Note  this  is  not  the  d  A  /<i  In  x  as  in  Equation  (l). 

The  lognormal  distribution  has  been  treated  at  leng'th  by  several 
authors  (AI-A3)  besides  Aitcheson  and  brown,  but  it  will  be  helpful 
to  summarize  some  of  its  essential  p'v.perties  here.  A  random  variable 


is  distributed  lognormally  if  the  logarithm  of  the  variable  is  dis¬ 
tributed  normally.  Such  a  variable  can  only  have  positive  values. 

The  value  of  the  variable  x  for  which  the  frequency  is  a  maixin.um  is 
the  modal  value  x  and  its  natural  logarithm 

In  X  = 

is  the  mean  value  of  the  In  x.  V/e  will  call  the  values  of  x  for 

which  t  =  +  1  by  x^^^  and  x  ^  At  these  values,  y  falls  to  exp  (-^)  of 

!L  2  2 

its  maximum  value  — - -  ,  -d  In  y/dt  is  unity,  and  d  y/dt  =  0. 

VSrt 

The  standard  deviation  a  is  given  by  In  x^^^  -  In  x  =  In  x  -  In  x 
the  difference  between  the  values  of  In  x  at  the  axis  of  symmetry  and 
the  points  of  inflection.  Thus  t  represents  the  departure,  in  o  units, 
of  In  X  from  its  mean  value. 

Finally, 

p  (x)  dx  =  -- —  exp 

V  2n 

=  P  (x,  X  +  dx) 

=  d  A(n.a^) 
d  In  X 

is  also  the  probability  that  In  x  lies  between  In  x  and  In  x  +  d  In  x. 

Figure  7  illustrates  several  presentations  of  the  same  distribution 
and  their  relation  to  one  another.  T^igure  "JA  is  a  logarithmic-probability 
presentation  for  Y  when  x  =  1  and  =  l/x  =  2.  Ibe  linear-probability 


1  1  In  x-d  1 

d  In  X 

"20 

0 

L  '  '  J 

presentation  of  the  normal  distribv/  obtained  by  letting  x' 


In  X 


is  shpwi  in  Figure  TB.  Here,  x  =  0,  x. ,  =  -x  ,  =  o  =  0.693'  This  is 

*“  ■**i  -  u 

converted  to  the  more  familiar  frequency  function  in  Figure  ^C  and 
to  the  less  familiar  Gaussian  plot  (a4).  These  in  turn  are  converted 
back  to  the  logarithmic  fom  in  Figure  70- 
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